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Abstract 

The study investigates the piezoelectric harvesting 
efficiency of ZnO microstructures with various structural 
configurations for electromechanical nanogenerators 
using finite element method simulation. The effects of 
ZnO structural parameters for a given cross-section 
area are simulated, including length, height, array 
number and shape. In addition, the harmonic 
responses and model frequencies are evaluated to 
analyze the natural mode shapes and force vibrations. 

The results show that the energy harvesting efficiency 
increases with decreasing ratio (Rh) of the diameter to 
the height. However, the trend of the harvesting 
voltage becomes almost saturated when the Raratio is 
less than 0.5.A single bulk structure of piezoelectric 
ZnO produces less micro-energy than that produced 
by an array structure for given cross-section area. The 
harvesting voltage increases with increasing number 
of arrays. The conversion efficiencies of 6, 11, 19 and 
22 arrays increases 93.4, 187.3, 345.6 and 537.2 %, 
respectively, that are compared to the single-bulk 
output powers. However, the energy harvesting 
efficiency deceases when the array number is over 33 
for arrays with a diameter of 14 \im. In addition, 
piezoelectric harvesting with ZnO arrays with a zigzag 
layer has higher efficiency than bulk structures. The 
results provide useful information for designing and 
fabricating micro-harvesting devices. 
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Introduction 

With increasing use of green energy, there is growing 
demand for higher efficiency energy harvestings 
[Yasar, 2012]. In general, solar energy, thermal energy, 
and mechanical vibrations [Philipps, Peharz, Hoheisel, 



and Hornung, 2010], [Felix, Ashcon, and Laurent, 
2012], [Xianzhi, Yumei, Ping, Yang, and Ming, 2011] 
can be harvested. Solar cells have been limited in 
low-light conditions despite their excellent power 
density. The thermal energy for harvesting is less 
efficient due to high power consumption. Mechanical 
vibrations can be converted into electrical energy 
using electromagnetic, electrostatic and piezoelectric 
devices. Generally speaking, piezoelectric devices 
have a high energy conversion for micro-harvesting 
[Roundy, Leland, Baker, Carleton, and Reilly, 2005] 
and have great promise for powering electronics and 
achieving self-powered electronic devices [Chen, Caofeng, 
Ying, and Zhong, 2012]. The piezoelectric materials 
such as zirconium titanate (PZT), aluminum nitride 
(A1N), barium titanate (BaTiOs), polyvinylidene 
fluoride (PVDF) and zinc oxide (ZnO) have been used 
for mirco-harvesting [McDonald, Genova, Rosenberg, 
Keating, Benedixen, and Junru, 2012], [Kumar and 
Kim, 2012], [Bai, Xu, Gu, Ma, Qin, and Wang, 2012]. 

In recent years, Al-doped ZnO nanorods have mostly 
attracted increasing attention due to their conductivity 
and transparency for applications in optoelectronic 
devices [Varadhaseshan and Sundar, 2012] [Wu, Wang, 
and Chen, 2007]. In addition, ZnO nanorods have an 
extremely stable field emission property [Xue, Li, Yu, 
Chen, and Yang, 2006] and have been considered for 
transparent conductive electrodes [Kong, Choi, Cho, Kim, 
Baek, and Lee, 2010] due to their high transparency, 
stability and high conductivity. According to reports, 
the harvesting energy of ZnO nanorods with 
Al-dopedis significantly gets enhanced at peak 
intensity of c-axis (002) orientation compared to the 
peak of hexagonal ZnO [Mondal, Kanta, and Mitra, 
2008]. 

In general, the piezoelectric harvesters are bulk [Bau, 
Ferrari, Tonoli, and Ferrari, 2011] or film [Jing, Miao, 
Xu, Olfatnia, and Norford, 2010] structures. The 
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piezoelectric materials exhibit an energy coupling 
between mechanical deformation and electrical energy 
when a normal stress is applied. Piezoelectric generators 
can utilize a cantilever beam, single nanorod or zigzag 
configuration to harvesting energy [Kanno, Ichida, 
Adachi, Kotera, and Shibata, 2012]. These structural 
configurations have a large aspect ratio and critical 
strain, high flexibility, and sensitivity to small forces. 
However, few studies have made simulation to 
determine the piezoelectric harvesting efficiency 
[Wang, Song, Liu, and Wang, 2007], [Xu, Qin, Xu, Wei, 
Yang, and Wang, 2010], [Zhu, Yang, Wang, and Wang, 
2010]. 

This study simulates piezoelectric harvesting efficiency 
of ZnO microstructural arrays with various structural 
configurations using finite element method (FEM) 
simulation. The effects of various ZnO structural 
parameters on a given cross-section areas are simulated, 
including length, array number and shape of ZnO. In 
addition, a modal analysis is conducted and the 
harmonic response of ZnO structures is evaluated to 
illustrate the excitation frequency at maximum 
deflection. 

Simulation Setup 

Harvesting Models 

Coventor Ware software (Conventor, Inc.) was used to 
simulate the harvesting efficiency of piezoelectric ZnO 
structures. A pair and an array of piezoelectric 
structures with a zigzag layer were designed to 
enhance bending and compression deformation, as 
shown in FIG. 1(a) and 1(b), respectively. The zigzag 
structure provides a normal stress and a bending 
moment to the array of piezoelectric ZnO, which enlarges 
the electric field due to higher electromechanical 
coupling. During the simulation, a half symmetry 
model of the ZnO structures was simulated to reduce 
memory use and time. The mesh elements of the 
rectangle and the cylinder were meshed as manhattan 
and tetrahedrous types, respectively. The mesh size, 
from 4 to 10 um, was in accordance with the size of the 
piezoelectric structures. For the boundary condition, 
the bottom of the ZnO substrate was considered as 
rigid and there was no deflection during simulation. 
A5umthickness rigid cover layer on the top of the 
piezoelectric structure was used to apply a load. The 
ZnO structures were assumed to be linear and to have 
a small deflection under an external load. The mass, 
damping and stiffness matrices were constant with 
time during the simulation. 



Zigzag Pressure 
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(a) A pair 

Pressure 




substrate 



(b) ZnO array rods with zigzag layers that enhance bending and 
compression deformation 

FIG. 1 STRUCTURES USED FOR PIEZOELECTRIC SIMULATION 

For modal and harmonic simulations of energy 
harvesting, rectangular, cylindrical, and 2x2, 4x4, 5x5, 
and 6x6 arrays of symmetric ZnO models were 
simulated, as shown in FIG. 2. All structures of the 
ZnO models hada length of 60 um and a total 
cross-section area of about 707 um 2 . The width of the 
rectangular structure and the diameter of the 
cylindrical structurewere 26.59 um and 30 um, 
respectively. The diameters of the 2x2, 4x4, 5x5, and 
6x6 arrays were 15, 7.5, 6and 5 um, respectively. The 
length, width and thickness of the substrate were 50, 
50, and 30 um, respectively, and the thickness of the 
top layer (not shown) was 5 um. 




FIG. 2 CHARACTERISTIC SIMULATION OF ZnO MODELS FOR 
PIEZOELECTRIC HARVESTING FOR A GIVEN TOTAL 
CROSS-SECTION AREA 

Piezoelectric Harvesting Mechanisms 

The mechanisms of piezoelectric harvesting are created 
using the electromechanical coupling characteristics. 
The mechanical properties, namely the mechanical 
stress (a) and mechanical strain (S), and dielectric 
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properties, namely the electric field (E) and the electric 
displacement (D), of the piezoelectric material were 
considered. The constitutive equation to describle the 
piezoelectric properties is based on d-type: 

s = s E cr + d t E 

\D = d G + e a E 



(1) 



where 
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s E , d, and e° are the compliance coefficient (m 2 /N), 
stiffness matrix (m/V) and dielectric constant (F/m), 
respectively, dt is the transpose of d. The superscript a 
indicates the dielectric constant at a constant mechanical 
stress. The superscript E indicates the compliance 
coefficient at a constant electric field. Assume that the 
ZnO crystal material is isotropic, orthotropic and 
polarized along the z direction. Therefore, the 
compliance coefficient, stiffness matrix and dielectric 
constant can be respectively expressed as: 
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Substituting equations (2) to (4) into (1) and representing 
the result in array form yields: 
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During the simulation, the 
were Sn=2.097xl0 5 , S12 = S 2 i= 



compliance coefficients 
1.211xl0 5 , S 22 =2.097xl0 5 , 



Sl3=S31=1.051xl0 5 , S23=S32=1.051xl0 5 , S33=2.109xl0 5 , S44= 

4.430xl0 4 and S55=S 66 =4.240xl0 4 MPa. The stiffness 
matrices were d25=-1.134xl0 5 , d26=-1.134xl0- 5 , d3i=-5.43x 
10- 6 , d32=-5.43xl0- 6 and d 2 4=0, d33=1.167xl0" 5 pC/(V-^m). 
The vacuum permittivity, density, and dielectric 
conductivity were £o=8.85xl0 12 C/(V-m), 5.68 g/cm 3 
and 2.5X10" 1 ps/|am, respectively. 

Results and Discussion 

The effects of ZnO geometric shape, a modal analysis 
and a harmonic response evaluation are illustrated in 
this section. 

Effect of Geometric Shape 

Most piezoelectric structures for energy harvesting are 
designed as rectangular cuboids due to the fabrication 
process. Rectangular and cylindrical structures with a 
cross-section area of about 707 |am 2 and a length of 60 
um were simulated at various shear stresses. The width 
of the rectangular and the diameter of the cylindrical 
were 26.59 \xm and 30 \xn\, respectively. The length, 
width and thickness of the substrate were 50, 50 and 
30 |am, respectively, and the thickness of top layer was 
5 \xm. The results show that the harvesting voltages 
and structural deflections of the cylindrical shape are 
higher than those of the rectangular cuboid for all shear 
stresses, as shown in FIG. 3. This is attributed to the 
fact that the cylindrical shape has a smaller moment of 
inertia than that of the rectangular cuboid. The 
moment of inertia means that an object resists 
rotational moment and torque about an axis. 
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FIG. 3 ENERGY HARVESTING CHARACTERISTICS OF 
RECTANGULAR AND CYLINDRICAL ZnO STRUCTURES FOR A 
GIVEN CROSS-SECTION AREA AND LENGTH 

FIG. 4 shows voltage versus shear stress for cylindrical 
structures with various diameters and a constant height 
of the ZnO cylinder for shear stresses of 1, 0.3, 0.5, 0.7, 
and 1.0 MPa. A constant height of the ZnO stucture is 
100 |am and the diameters of the ZnO stucture are 20, 
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40, 60, 80, 100 and 120 [im. The geometric ratios, Rh= 
cf)/ h, of the diameters, cf), and a constant height, h, are 
0.4, 0.6, 0.8, 1.0, and 1.2 for energy harvesting. The 
harvesting voltage increases with decreasing 
geometric ratio Rh. The harvesting voltage versus ZnO 
diameter for a constant pressure of 1 MPa is shown in 
the inset of FIG. 4. The harvesting voltage curve 
becomes almost saturated when the diameter of the 
piezoelectric ZnO is larger than the height of structure. 
The curve of the harvesting voltage is simillar to 
negative exponential function. 
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FIG. 4 VOLTAGE VERSUS SHEAR STRESS FOR CYLINDRICAL 
ZNO STRUCTURES WITH VARIOUS DIAMETERS AND A 
CONSTANT HEIGHT OF ZnO CYLINDER FOR VARIOUS SHEAR 
STRESSES 
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FIG. 5 VOLTAGE VERSUS SHEAR STRESS FOR CYLINDRICAL 
STRUCTURES WITH VARIOUS HEIGHTS AND A CONSTANT 
DIAMETER OF ZnO CYLINDER FOR VARIOUS SHEAR STRESSES 

FIG. 5 shows voltage versus shear stress for cylindrical 
structures with various heights and a constant 
diameter for shear stresses of 1, 0.3, 0.5, 0.7 and 1.0 
MPa. The ratios Rd is defined as the heights (h) of ZnO 
cylinder anda constant diameter (cf)) of ZnO cylinder. 
The ratios, Rd=h/cf), of ZnO cylinder for energy 
harvesting simulation are 2, 1.67, 1.0, 0.67, 0.5 and 0.33. 
The diameter of the ZnO stucture is 60 um and the 
heights of the ZnO stuctureare 20, 30, 40, 60, 100 and 



120 um. The results show that the harvesting voltages 
of the ZnO stucture increase with increasing geometric 
ratioRd. The harvesting voltages for R equal 2, 1.67, 1.0, 
0.67, 0.5 and 0.33 at a shear stress of IMPa are 6.9, 5.23, 
2.2, 1.27, 0.49 and 0.48 mV, respectively. The conversion 
efficiencies of Rd= 0.33, 0.5 and 0.67 decreased 78.2, 
77.7 and 42.3 %, respectively, compared to the 
geometric effect of Rd=l. However, the conversion 
efficiencies of Rd=1.67 and 2.0 increased 137.7 and 
213.6%, respectively, also compared to the geometric 
effect of Rd =1. The harvesting voltage trend of the 
ZnO becomes almost saturated when the Rd ratio is 
less than 0.5. It is can be predicted that there has small 
electric diploemoment when the Rd is less than 0.5. 

Modal and Harmonic Analysis 

In general, the harvesting mechanism of a piezoelectric 
structure is produced from structural resonances due 
to vibrations at the excitation frequency. An modal 
analysis can be used to determine the natural mode 
shapes and the resonant frequencies of the ZnO 
structures for optimal energy harvesting during free 
vibration. FIG. 6 shows the ZnO modal analysis results 
of a rectangular, a cylindrical, and 2x2, 4^4, 5x5, and 
6x6 array structures with the same total cross-section 
areas. Simulation results show that the rectangular 
ZnO structure has a lower resonant frequency than 
that of the cylindrical structure at different modal 
types. It is supposed that a rectangular has larger 
second moment of area than a cylinder. A ZnO 
structure with fewer arrays has a higher resonant 
frequency. Therefore, a structure with a low resonant 
frequency is more efficient for energy harvesting due 
to a less amount of driving energy. The energy 
harvesting for microstructures at high-frequency MHz 
range is not suitable for the use of array structures. 




Modal type 



FIG. 6 ZnO MODAL ANALYSIS RESULTS FOR A RECTANGULAR, 
A CYLINDRICAL, AND 2x2, 4x4, 5x5 AND 6x6 ARRAY 
STRUCTURES FOR A GIVEN TOTAL CROSS-SECTION AREA 
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The time-dependent equations of ZnO motion 
undergoing steady-state vibration were analyzed. All 
loads and deflections sinusoidally vary at a given 
frequency. The harmonic analyses can be used to 
obtain the deflection amplitudes of ZnO vibration at 
given points in the structures as a function of forcing 
frequency. FIG. 7 shows the harmonic deflections of 
various ZnO models for various harmonic forces at a 
constant resonant frequency. The ZnO model with a 
6><6 array had the highest deflection for a given 
harmonic force. Energy harvesting efficiency of ZnO 
array increases with the increasing deflection. The 
deflections of a rectangular, a cylindrical, and 2^2, 4x4, 
5x5 and 6x6 array structures at a harmonic force of 500 
uN were 0.708, 0.943, 4.215, 6.759, 9.786 and 14.748 urn, 
respectively. Therefore, ZnO structures for energy 
harvesting should be designed as an array of cylinders. 
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FIG. 7 HARMONIC DEFLECTIONS OF ZnO MODELS FOR A 
GIVEN CROSS-SECTION AREA AND VARIOUS HARMONIC 
FORCES AT A CONSTANT FREQUENCY 

Array Analysis 

The energy harvesting of single bulk and array 
structures with a given cross-section area and a height 
of 100 um were investigated under various shear 
stresses. The diameter and the height of the single bulk 
were 60 and 100 um, respectively. The diameters of 
ZnO with 6, 11, 19, 22 and 33 arrays were 30, 24, 19, 16 
and 14 um, respectively. The harvesting voltages of a 
single bulk are smaller than those of the arrays, as 
shown in FIG. 8. In addition, the harvesting voltage of 
ZnO array increased with increasing number of arrays. 
The voltages for 1, 6, 11, 19 and 22 arrays are 0.543, 
1.05, 1.56, 2.42 and 3.46 mV, respectively. The 
conversion efficiencies of 6, 11, 19 and 22 arrays 
increase 93.4, 187.3, 345.6 and 537.2 %, respectively, 
compared to the single bulk output powers. The 
harvesting energy of ZnO decreased when the number 
of arrays was over 33. The structure rigidity increases 
with increasing number of arrays. 
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FIG. 8 PIEZOELECTRIC HARVESTING RESULTS FOR A SINGLE 
STRUCTURE AND 6, 11, 19, 22, AND 33 ARRAYS 

Conclusion 

The energy harvesting characteristics of zinc oxide for 
microharvesters were investigated. The harvesting 
voltage curve becomes almost saturated when the 
diameter of the piezoelectric ZnO is larger than the 
height of the structure (about 100 um). A rectangular 
ZnO structure has a lower resonant frequency than a 
cylindrical structure for a given cross-section area at 
different modal types. The harvesting voltage of ZnO 
increased with increasing number of arrays. However, 
the harvesting energy of ZnO decreased when the 
number of arrays was over 33. The results show that ZnO 
can be applied to electromechanical nanogenerators. 
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